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ABSTRACT 

We study the relationships between galaxy environments and galaxy properties related to disk 
(re)growth, considering two highly complete samples that are approximately baryonic mass limited 
into the high-mass dwarf galaxy regime, the Environmental COntext (ECO) catalog (data release 
herein) and the B-semester region of the REsolved Spectroscopy Of a Local VolumE (RESOLVE) sur¬ 
vey. We quantify galaxy environments using both group identification and smoothed galaxy density 
field methods. We use by-eye and quantitative morphological classifications plus atomic gas content 
measurements and estimates. We find that blue early-type (E/S0) galaxies, gas-dominated galax¬ 
ies, and UV-bright disk host galaxies all become distinctly more common below group halo mass 
~ 10 11 ' 5 Mq, implying that this low group halo mass regime may be a preferred regime for significant 
disk growth activity. We also find that blue early-type and blue late-type galaxies inhabit environ¬ 
ments of similar group halo mass at fixed baryonic mass, consistent with a scenario in which blue 
early types can regrow late-type disks. In fact, we find that the only significant difference in the 
typical group halo mass inhabited by different galaxy classes is for satellite galaxies with different 
colors, where at fixed baryonic mass red early and late types have higher typical group halo masses 
than blue early and late types. More generally, we argue that the traditional morphology-environment 
relation (i.e., that denser environments tend to have more early types) can be largely attributed to 
the morphology-galaxy mass relation for centrals and the color-environment relation for satellites. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution - ultraviolet: galaxies 


1. INTRODUCTION 


For decades, astronomers have observed that the prop¬ 
erties of galaxies in the local universe, including ap¬ 
pearance, star formation rate, and gas content, depend 
on the surrounding environment (e.g., as rev iewed by 
Bower fc Balogh||2004 iBoselli fe Gavazzl||2006 ). G 
ies of different morphological types, in particular, 
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Dressier (JIDbU) reported the so-called 


increase with increasing environmental densit y wit hin 
rich clusters wh ile spiral fractions decrease. |Postman 


& Geller (1984) showed that a similar relationship be- 


tween galaxy morphology and local density also holds 
in the lower-density group environment. However, sev¬ 
eral authors found the conflicting result that significant 
variation in morphology with environment exists only in 
the richest clusters (e.g., |Maia fe da Costa 1990 Whit¬ 
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vironments has since been corrobor ated by a variety of 


other authors (e.f 


Tran et al. 2001 

Helsdon & Ponman 

2003 Goto et al. 
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03} Calvi et al. 2012). 


galaxies considered can have a significant impact on the 
observed form of the morphology-density relation (e.g. 


alvi et alT] 20121 IWilman fe Erwin| 

20121). While brinkwater et al. (|200i|) observe a tradi- 

tional morphology-density relation m the dwarf galaxy 
population within a single cluster, other studies find 
that the morphology-density relation takes on different 
forms for lo w-to-intermediate mass galaxies. |Kannappan| 
et al. (2009 hereafter KGB) find suggestive evidence that 
low-to-intermediate stellar mass E/SOs occupy low den¬ 
sity environm ents s imilar t o tho se of spirals at the same 
masses, while Calvi et al. (2012) find that morphologies 
for intermediate mass galaxies are not closel y related to 
environment, except in clusters. Interestingly, |Hogg et al.] 
(2003) find that the mean environmental overdensity for 
red, typically assumed to be early-type (E/S0), galaxies 
reaches a minimum at intermediate masses/luminosities 
(around L*). The conflation of morphology and color, 
of course, complicates interpretation of the Hogg et al. 
result, s ince both optically red, “passive” spiral galax- 
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Dressier 


galaxies (e.g., KGB; Schawinski et al. 2009|) are known 
components of the galaxy population. 

It is likely that the existence of such non-traditional 
color-morphology pairings in the galaxy population is 
an important factor driving the observation that trends 
in galaxy color (or more direct star formation property 
measures) versus environment can differ from observed 
morphology-density trends. Beginning with the analy- 
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Fig. 1.— The ECO catalog region in sky c oord inates, with the ECO mass-limited sample represented by points color-coded according to 
group halo m asses estimated as d escribed in §3.6| The RESOLVE-A region is outlined in red, and the region of overlap with the ALFALFA 
a.40 catalog (Haynes et al. 201 lj is indicated by the purple crosshatched strips. 


sis of Kennicutt (1983), it has often been observed that 
cluster galaxies show typically lower levels of star for¬ 
mation than galaxies in less rich environments. Several 
authors have subsequently found that such star forma¬ 
tion or color trends with environment cannot be com- 


trends (e.g., Koopmann & Kenney 

1998; Lewis et al. 

2002 

Ghristfein & Zlabiudoli 2005; V> 

/eiikaia et af. 2009). 


Others have taken this conclusion further, interring that 
trends in star formation/color are actually more closely 
linked to environmental conditions than are trends in 


structural/morphological parameters (e.g., 

Blanton et al. 

2003; Kauffmann et al. 

2004; 

Blanton et a 

. 200m Quin- 

tero et al. 2006 Bamford et a 

. 2009[ Skibba et af. 2009). 


observed to share a close link to both star formation and 
the ambient envi ronmental condit i ons a round galaxies. 


First observed by Davies & Lewis (1973) for galaxies in 
the Virgo Cluster, the result is now well established that 
cluster galaxy populations typically display lower levels 
of atomic gas than do sim ilar popula tions in le ss d ense 


environments (e.g., Giovanelli fc Haynes 


___ riaynei 

et al.||1984( |Giovanelli fc Haynes||1985| Havazzi|| 1987[ 

Bravo-Alfaro et al.||2000| ISolanes et al.|2001| see also re- 


1983 


Pay: 
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view by van Gorliomj 


are typically observec 


n||20C)4j). 

a to Eav 


Likewise, cluster galaxies 
have HI gas disks that are less 


extended than those in lower density environments (e.g., 


Giovanelli & Haynesl 1983 

Warmels 1988| Cayatte et al. 

1994 Bravo-Alfaro et al. 

2000p. Galaxy gas properties 


lower density environments as w ell, e.g., a s shown in the 


simulations of Kawata & Mulchaey (2008) where “stran¬ 
gulation” or stripping of a hot galaxy halo gas compo¬ 
nent, which could otherwise cool to provide a cold gas 
reservoir, in a larger potential is effective in a group envi¬ 
ronment with halo mass ~ 10 12 ' 9 Mq. Still further down 
the environmental density scale, gas-rich galaxies have 
been observed to be one of the most weakly clustered 
galaxy populations, th at is, typically found in the lowest 
density environments (Basilakos et al.|2007 Meyer et al. 


2007 Martin et al. 

2012 

Li et al. 

2012) 


retical studies that gas accretion into galaxies, whether 
in a “cold” or “hot” mod e, is most effective where group 
halo masses are low (e.g., Birnboim fc Dekel|2003 Keres 
et al. H20051 IDekel &; Birnboim||2006| |Keres et al. Pl)U9| 


Nelson et al.||2013[). 


In the prevailing hierarchical model of galaxy evolu¬ 
tion, galaxies are thought to experience morphological 
transformations not just from late to early type, for ex¬ 
ample through merging/quenching processes, but also 
potentially from early to late type, through a disk re- 
qrowt h process that may be enabled by gas accretion 


.g., |B arnes 2002; |Stei 
nato et al.||2007p. Opservationally, the operation of such 


Steinmetz & Navarro 


2002 Gover- 


a disk regrowth process is difficult to confirm. 

However, several hints have recently emerged that such 
a scenario is plausible. One such hint lies in the exis¬ 
tence of the blue or “blue-sequence” E/S0 population, 
consisting of morphologically early type galaxies that lie 
on the blue sequence in color-stellar mass space. Blue- 
sequence E/SOs are typically found in non-cluster en¬ 
vironments and exist primarily at stellar masses less 
than ~ 10 10 5 Mq (the bimodality mass of Kauffmann 


et al.| 2007 hereafter Mb) but are most common below 
the “ga s-richness threshold ” stell ar mass of ~ 10 9 7 Mq 
(KGB). Kannappan et al. ((2013 hereafter K13) argue 
that the bimodahty and gas-richness threshold mass 
scales mark two distinct transition points between galaxy 
refueling regimes, with galaxies below the threshold scale 
typically experiencing high levels of external gas accre¬ 
tion and stellar mass growth. Consistent with this pic¬ 
ture, low-mass blue E/SOs contain sufficient gas reser¬ 
voirs and specific star formation rates to allow the growth 
of evolutionarily significant disk struct ures on relatively 
short timescales (KGB; Wei et al.|2010). Another hint of 
disk regrowth in E/SOs is the observation of extended UV 
emission, associated with recent star formation, around 
a number of nearby galaxies with early-type morphol- 
ogy (e.g., Donovan et al. ||2009 Cortese & Hughes 2009 
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Fig. 2.— The ECO catalog region in RA vs. line-of-sight distance coordinates, in slices of ~5 degrees in Dec, increasing from the bottom 
left panel. Galaxies included in the final ECO mass-limited sample analyzed here are indicated by dots color-coded according to their 
group halo masses as in Fig.fl] while galaxies outside this sample but present in our merged parent redshift catalog are indicated by grey 
dots. The outer limits of the JECO catalog “buffer” region are outlined in red, and the region of ECO we consider interior to this buffer is 
indicated in green. Black lines in the upper left and bottom right panels indicate the approximate line-of-sight extent of two groups that 
extend significantly outside the ECO region and are subject to boundary completeness correction factors as described in j |3.8| 


Thilke r et al. 2010 Moffett et al. 2012). Moffett et al. 
j2012|) observe that such extended UV structures are 
common in low-to-intermediate mass early-type galaxies 
and that a particular class of “UV-Bright” or UV-B disk 
galaxies is marked by a high potential for ongoing star 
formation. Linking these two populations together, UV- 
B disks are also strongly associated with the low-mass, 
blue-sequence E/SO population, supporting the ide a that 
these galaxies may be engaged in disk regrowth. |Stark 


et al. (2013) report evidence for gas as well as stellar disk 


regrowth in low-mass, post-starburst E/SO galaxies. 

In this contribution, we employ two unusually com¬ 
plete volume-limited galaxy samples, both of which ex¬ 
tend into the “high-mass” dwarf galaxy regime (reaching 
baryonic masses ~ 1O 9 - 3 M 0 ), to probe disk (re)growth in 
a variety of environments. We seek to answer three major 
questions. (1) Does environment play a role in enabling 
gas or stellar disk growth? (2) Does the morphology- 
density, or more generally morphology-environment, re¬ 
lation behave as might be expected if disk regrowth is 
effective in transforming galaxy morphology? In particu¬ 
lar, are the typical environments for blue early- and late- 
type galaxies similar in the galaxy mass regimes in which 
disk regrowth occurs? (3) Are there specific group halo 
mass scales implicated in evidence for disk regrowth? 

We address these questions in part by examining the 
detailed form of the morphology-environment relation, 


including possible variations with galaxy mass scale 
and central/satellite designations. In addition, we con¬ 
sider an alternative way of formulating a morphology- 
environment relation. If the traditional formulation of 
the morphology-environment relation can be considered 
to quantify the probability of a galaxy exhibiting a par¬ 
ticular morphology given some environment, P(M|E), 
then an alternative way to frame this relation is to quan¬ 
tify the probability of a galaxy inhabiting a particular 
environment given its morphology, P(E|M). This alter¬ 
native formulation provides a useful framework for un¬ 
derstanding the typical environments of galaxies with 
different morphologies. We further examine the typi¬ 
cal environments of different classes of galaxies linked to 
disk growth, including blue-sequence early types, galax¬ 
ies with substantial atomic gas reservoirs, and early-type 
galaxies that display recent UV-detected disk star for¬ 
mation. 

We begin by introducing the galaxy samples under 
consideration in f[2] In M we describe our main data 
analysis methods. In ^4 we report a variety of results 
from this analysis, including the observed forms of both 
the traditional and alternative morphology-environment 
relations, in particular finding that the morphology- 
environment relation disappears for low baryonic mass 
central galaxies. We also find that blue-sequence early 
type, gas-dominated, and disk-growing populations rise 
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Fig. 3. — Illustration of the selection of the ECO catalog. Panel a shows the M r and Mf, ary distributions of the initial ECO catalog (dots 
and purple inset Mf, ary histogram), where the horizontal red line indicates the M r < —17.33 redshift completeness limit and the vertical 
green lines indicate the final mass cut we adopt to create an approximately baryonic mass limited sample to M> }ary > 10 9 ' 3 A/q. Panel 
b shows the magnitude and cz limits of ECO, where light grey dots indicate the M r < —17.33 sample and black dots indicate our final 
approximately mass-limited sample with M^ ary > Ilf' 5 Mq and group membership within the ECO volume. 

not present in either one to be recovered through in¬ 
clusion of the other, with SDSS typically providing red- 
shifts for fainter objects than the UZC. Though defined 
by the overlap of these two catalogs, the ECO cata¬ 
log also incorporates redshifts from the REsolved Spec¬ 
troscopy Of a Local Volum E (RE SO LVE; Kannappan et 
, HyperLEDA ([Paturel et al.[ 


350 


300 


Mass lim. ECO 
0 ECO with SDSS redshift 



al. 


m prep. 


Driver et al. 2011), 2dF (Coiless et al. 


.'2003), GAMA 
2001 ). and 6dF 


(Jones et al.|2U09) surveys. The ECO region was also se- 


Fig. 4. — Illustration of the additional completeness of the fi¬ 
nal Mb ary -limited ECO catalog over the SDSS redshift catalog. 
The M r distributio n of the ECO sample in our reprocessed mag¬ 
nitude system (see j ]3.1| for details) is indicated by the purple solid 
histogram, and the distribution of ECO galaxies that have SDSS 
redshifts is indicated by the black hashed histogram. 

in prominence in environments with M^aio < 10 11 ' 5 Mq. 
In we show that the forms of both the traditional 
and alternative morphology-environment relations we 
observe are consistent with expectations of the disk re¬ 
growth model and discuss the idea that the low group 
halo mass regime below ~ 10 115 Mq appears to be a 
preferred regime for disk growth. Finally, we provide a 
brief summary of our major results in f|6] 


2 . 1 . 


2. SAMPLES 
ECO catalog 

The ECO, or Environmental COntext, catalog is the 
largest sample we consider and includes the greatest di¬ 
versity of galaxy environments (see Table 1 for detailed 
ECO galaxy properties). The ECO catalog region was 
chosen as the largest contiguous region on-sky whe re the 
highly complete Updated Zwicky Catalog (UZC; Falco 


lected to enclose the RESOLVE A-semester survey vol¬ 
ume plus a minimum 1 Mpc “buffer” in all directions 
(see sky coverage in Fig. [lj). This buffer region, cho¬ 
sen with a size comparable to typical group halo virial 
radii at the present epoch, exists to mitigate potential 
edge effects in calculating galaxy environment metrics, 
such that only galaxies in the buffer region should have 
environmental measures strongly affected by the loss of 
nearby galaxies outside the catalog boundaries. The far 
side limit of 7470 km/s was selected to encompass both 
the aforementioned 1 Mpc (equivalent to 70 km/s for 
Ho = 70 km s -1 Mpc -1 ) buffer beyond the RESOLVE 
cz limit of 7000 km/s and an additional allowance to 
compensate for group peculiar velocities up to 400 km/s. 
The near side buffer zone cz limit of 2530 km/s was simi¬ 
larly chosen to expand the ECO volume as much as pos¬ 
sible beyond the near-side RESOLVE cz limit of 4500 
km/s while avoiding the effects of Virgo Cluster region 
velocity-space distortions (see Fig.[2|. We consider the 
velocity limits of the non-buffer ECO volume to be 470 
km/s away from the buffer edges, or 3000 km/s < cz < 
7000 km/s. 

The ECO catalog represents a cross match between 
sources with measured redshifts found in the UZC, SDSS 


(including data releases 6, 7, and 8 ; lAdclman-McCarth; 


et al. 2008 


Aihara et al. 

and bdf* 


2011 


f 


cata- 


et al. 


et al. 


1999) and Sloan Digital Sky Survey (SDSS; York 


20001 redshift databases overlap, allowing objects 


Abazajian et al. 2009 
HyperLEDA, RESOLVE, GAMA, 2HF 
logs with a 15” matching radius on sky. New sources 
are added to ECO from each of the constituent cata¬ 
logs whenever they do not match to a previously in¬ 
cluded ECO source. The resulting catalog has also been 
inspected by eye for duplicate entries caused either by 
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“shredding” of SDSS photo metric ob jects into multiple 
galaxy pieces (as described in Abazajian et al. 2004) or by 
centering/coordinate errors occasionally larger than the 
cross-matching radius. Such duplicate entries, making 
up ~5% of the galaxies originally considered for inclu¬ 
sion, have been removed from our catalog. As illustrated 
in Fig. [4j the majority of ECO galaxies are present in the 
SDSS redshift survey, but the number of galaxies added 
from other sources is significant, at approximately 7% of 
the final catalog. 

Initially, sources with measured positions and redshifts 
inside the ECO volume are considered for potential mem¬ 
bership regardless of any previously measured catalog 
magnitudes they may or may not have. We then use 
custom photometric measurements performed on SDSS 
i mag ing f rames for all such potential ECO members (see 
(3.1 and Eckert et al. 2015) to determine a defining 
magnitude limit for the ECO catalog. The complete¬ 
ness limit of the SDSS redshift survey at 7000 km/s is 
M r = —17.23 (in DR7 catalog Petrosian magnitudes cor¬ 
rected for foreground extinction), and our reprocessed 
magnitudes are typically brighter than the SPS S cata¬ 
log values by approximately 0.1 mag (see §3.1| for de¬ 
tails). Thus, a potential completeness limit for the ECO 
catalog motivated by the SDSS completeness would be 
M r < —17.33, which we use t o produce an initial group 
catalog as described in §3.6.1| However, we seek to cre¬ 
ate a baryonic mass limited final sample for our analysis, 
where baryonic mass (Mb ary ) is defined here as stellar 


plus ato mic gas mass and estimated as described in §3.2| 
and |3.7| As illustrated in Fig. [3ji, within the confines 
of the M r = —17.33 magnitude completeness limit, we 
can construct an approximately baryonic mass limited 
sample with Mt mr y > 10 9 3 Mq while leaving aside only 
a relatively small number of high mass-to-light ratio ob¬ 
jects. 

The final ECO sample we analyze meets both the above 
limit and the additional criterion that th e cen ter of the 
group to which each galaxy belongs (see (3.6 for details 
of group membership determination) must lie within the 
limits of the non-buffer region of ECO, that is, the group 
center must have 3000 km/s < cz < 7000 km/s and 
RA/Dec > IMpc from the buffer edges on sky at its 
redshift (see Figs. [ 2 ] and [3 }d)[^] This final sample contains 
6716 galaxies. 

Two partially overlapping subregions of the ECO cat¬ 
alog are given special attention in this paper. The 
“ECO+A” region is defined by the overlap of the ECO 
catalog and public Arecibo Legacy Fast ALFA (AL¬ 


FALFA) a.40 catalog ( [Haynes et al. 


resents the portion of the ECO samp 


2011) and rep- 
e for which di¬ 


rect HI mass determinations are available and used in 
our analysis (see Fig. [lj. The ECO+A region en¬ 
compasses the full range of environments found in the 
ECO sample. The “ECO+G” subsample is defined 


by the availability of archival GAL EX imaging (Mor¬ 


rissey et al. 2007 GALEX MAST GR6/7 archive at 


http://galex.stsci.edu/GR6/) with exposure times 
> 1000s in the NUV band, sufficient to detect extended 


7 Note that we do not recalculate galaxy M r values by consid¬ 
ering galaxies to lie at their group center redshifts. We find this 
recalculation would make a negligible difference in our overall sam¬ 
ple membership. 


ECO+G subsample a 
Full ECO 
RESOLVE-B B 



12 13 

Log( M hal0 [M sun ]) 

Fig. 5. — The distribution of galaxies by group halo mass in 
the mass-limited ECO, ECO+G, and RESOLVE-B sam ples, with 
completeness correction factors applied as described in §3.8| The 
red dashed vertical line indicates the group halo mass completeness 
limit for ECO, and the RESOLVE-B frequency histogram has been 
rescaled to match ECO at this limit. We have selected the ECO+G 
sample to have a similar environment distribution to the full ECO 
sample. 

UV disk structures (e.g., |Thilker et al.||2007|). Since 
GALEX imaging coverage over the EGO sky region is 
patchy, we select multiple fully covered subregions of 
ECO as our final ECO+G subsample, largely coincident 
with RESOLVE-A but extending to larger Dec and in¬ 
cluding several slices through rich clusters. Together 
these regions closely reflect the full ECO environment 
distribution (see Fig. [5]). 

2 . 2 . RESOLVE-B 

The B-semester region of the RESOLVE survey, which 
covers most of the SDSS “Stripe 82” region, is used as 
a comparison sample in this analysis. The RESOLVE- 
B region environment distribution is illustrated in Fig. 

t where primary differences compared to ECO are 
e lack of M halo > 10 13 5 M o groups while M halo ~ 
10 13 Mq groups are overrepresented in RESOLVE-B. The 
RESOLVE-B subsample has the advantage of greater 
completeness than the ECO catalog, due to deeper SDSS 
imaging and redshift coverage, plus the further redshift 
completion efforts of the RESOLVE survey. By compar¬ 
ison with this extra-complete sample, we assess the ef¬ 
fects of incompleteness in the ECO sample and derive 
completeness corrections that can be applied to ECO 
(median correction ~1%; see ( |3.8| ). For all galaxies in 
RESOLVE-B, morphological classification has been per¬ 
formed by a team of classifiers, providing both the basis 
for calibration of, and a comparison to, the quantita¬ 


tive morphological classifications used for ECO (see (3.4 


RESOLVE classifications details in Kannappan et al., m 
prep.). 


3. METHODS 

In this section, we describe our methods of custom 
photometric processing, galaxy color and stellar mass 
estimation, galaxy morphology and UV disk classifica¬ 
tion, environment metric calculation, atomic gas mass 
estimation, and correction for the incompleteness of the 
ECO sample. Throughout our analysis, we calculate 
distances according to D = cz/H 0 and take H 0 = 
70 km s _1 Mpc~ x unless otherwise noted. We estimate 
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Fig. 6. — Color vs. stellar mass for the ECO catalog sample, 
where (u — r) e represents an internal e xtinc tion corrected color 
derived from the SED fitting code (see j ]3.2| |. Dark grey points 
indicate individual ECO galaxies in our mass-limited sample, light 
pink points indicate ECO galaxies that do not enter this final sam¬ 
ple, and density contours of the mass-limited sample distribution 
are shown in purp le. T he green line indicates our chosen red/blue 
sequence divider ( j|3.3||. 


binomial confidence intervals on population proportions 
according to the Bayesian approach of Cameron (2011). 


3.1. Imaging/Photometry 

As a result of the problematic nature of obtaining ac¬ 
curate estimates of galaxy properties with bulk SDSS 
pipeline-processed data and our desire to study spatially 
resolved parameters not necessarily computed in catalog 
data products, we have undertaken a custom reprocess¬ 
ing of SDSS, 2MASS, and GALEX imaging for all ECO 
sample galaxies considered here. The ECO photomet¬ 
ric reprocessing mimi cs the metho ds de veloped for the 
RESOLVE survey (see Eckert et al.|2015 for full details). 

To summarize this reprocessing: we retrieve imaging 
frames in ugriz , JHK , and NUV b ands via a utomated 
queries to the SDSS D R8 (Aihara et ah 2011), 2MASS 
(Skrutskie et al. 2006|) , and GALEX (Morrissey et al. 


20071 archives respectively. Greater than 99% ol our 


galaxies are covered by 2MASS and ~30% by GALEX. 
Photometric processing proceeds first on the SD SS imag¬ 
ing, where SExtractor ( Bertin &; Arnouts||1996 ) is called 
to identify sources from an r-band image and create a 
corresponding mask image wherein sources other than 
the target are masked. SExtractor parameters for the 
target galaxy are then used as an initial input to the 
IRAF ellipse task, which fits the galaxy isophotes in a 
co-added gri image while allowing the PA and ellipticity 
to vary. From this free ellipse fit, an optimal ellipticity 
and PA corresponding to the outer disk is determined. 
A fixed ellipse fit is then performed, using these outer- 
disk parameters, on the images in each band individu¬ 
ally. Total magnitudes are determined from the resulting 
profiles by several methods: large aperture summation, 
exponential disk fitting, curve-of-growtlr extr apo lati on, 
and outer-disk color correction (see |Eckert et al. 2015). 
Comparing the results of these methods yields an esti¬ 
mate of the systematic errors, which are combined with 
the purely photometric errors to obtain the final magni¬ 
tude error estimates. 

The automated SExtractor masking procedure has 
been tuned to give reasonable results for the majority 


of galaxies, but it is possible for the automatically gen¬ 
erated masks to either mask parts of the galaxy under 
consideration or fail to mask nearby sources not asso¬ 
ciated with the target galaxy. To identify potentially 
problematic masks, we flag objects for which the mag¬ 
nitude estimation procedure has failed, the extracted r- 
band profile signal does not rise above eight times the 
sky noise, or the r-band profile does not extend beyond 
the calculated r-band 90% light radius. These mask im¬ 
ages flagged as potentially problematic are inspected by 
eye and edited by hand where necessary to better reflect 
distinctions between the target galaxy and other nearby 
sources. 


3.2. Color and Stellar Mass Estimation 

Using the full complement of total magnitudes, includ¬ 
ing NUV where available, galaxy stellar masses are esti¬ 
mated from a spectral energy distribution (SED) fitting 
procedure. We use a recently updated versi on of the 
likelihood-based stella r mass estimation code of |Kannap^| 
pan & Gawiser (2007), which is described fully by K13. 
This procedure uses a suite of composite stellar pop ula¬ 
tion models co nstructed from old and young Bruzual & 
Chariot (2003) stellar populations. These model stellar 
populations are combined in various fractions by mass, 
with 13 allowable young population fractions of 0.001, 
0.002, 0.005, 0.011, 0.025, 0.053, 0.112, 0.220, 0.387, 
0.585, 0.760, 0.876, and 0.941. Four possible metallic- 
ities (Z = 0.004, 0.008, 0.02, and 0.05) and 11 possible 
extinction values (tv = 0, 0.12, 0.24...1.2) are also used. 
The young population model grid is constructed to simu¬ 
late both continuous and bursty star formation histories 
by including models with constant star formation from 
1015 Myr in the past to various end points 0-195 Myr 
in the past and simple stellar populations (SSPs) with 
ages 360, 509, 641, 806, and 1015 Myr. The old popu¬ 
lation model grid includes SSPs with ages 2, 4 6, 8, lt^ 


and 12 Gyr. A Chabrier initial mass function (Chabrier 
2003) is used in these calculations, which yield a s teU 
far mass zer o point consistent with that of Kauffmann 
et al. (2003). The code also outputs internal extinction 
corrected model fit colors, which more cleanly separate 
the red- and blue-sequence galaxies than raw measured 
colors, and we designate these colors with a superscript 
“e”. We also use model fit colors without any internal 
extinction correction, denoted with superscript “model”. 


3.3. The Red and Blue Sequences 

To separate red- and blue-sequence galaxies, we choose 
a dividing line between the red and blue color stellar 
mass loci defined by our extinction-corrected colors and 
stellar mass estimates (see Fig. J6|. We determine this 
division based on double Gaussian fits to the red- and 
blue-sequence color distributions in two high and low 
stellar mass regimes where the sequences are well de¬ 
fined (log M* < 9.5 and log M* > 9.5). Our divider is 
then defined by the color halfway between the fit peaks 
and the median stellar mass in each mass regime. The 
slanted divider in the intermediate mass regime is defined 
by the line connecting these two Gaussian-fit-determined 
points. The equation of the divider is: 

( 1.457 logM* < 9.1 

(u-r) e = \ 0.24 x logM* - 0.7 9.1 < logM, < 10.1 (1) 

11.7 logM* > 10.1. 
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Fig. 7. — Quantitative morphology metrics applied to by-eye classified galaxies in the RESOLVE and ECO samples (see i |3.4| for details). 
Panel a shows a cut in C r vs. M r applied to ECO and RESOLVE galaxies classified by eye as early and late types. The solid green line 
shows an optimized morphology discriminant in this parameter space ( C r = —0.2 X M r — 1), which performs poorly in duplicating the 
by-eye morphological classes. Panel b shows the distribution of these same galaxies in the pa vs. M r parameter space. The solid green lines 
show our optimized morphology discriminant in this parameter space, which gives improved classification error rates over the concentration 
index approach. 


3.4. Morphology Classification 

To calibrate a quantitative morphology cut for appli¬ 
cation to the ECO galaxies, we use by-eye morpholog¬ 
ical classifications from the RESOLVE survey, the A- 
semester sample of which is largely a subset of the ECO 
catalog (see Kannappan et ah, in prep, for full details). 
RESOLVE galaxies that were given uncertain classifica¬ 
tions by the classifiers are omitted from consideration. 
Since this comparison sample has just over 1000 galax¬ 
ies, considerably fewer than the full ECO sample, and 
possesses few very bright galaxies, we also add to our 
morphology calibration sample those galaxies in ECO 
th at have been prev i ously classified by eye in the catalog 
of Nair fe Abraham (2010) or by the Galaxy Zoo Project 


(Lintott et al. 2011). We use only the “clean and de- 
biased” Gala xy Zoo classifications referenced by |Lintott| 
et al.| (20111, which require 80% of classifiers to agree 
on the chosen morphological type and debiasing with re¬ 
spect to the effects of lumino sity, size, and dista nce on 
the classifications (see §3.1 of Lintott et al., 2011). 

Based on comparisons to this by-eye classified sam¬ 
ple, an optimal quantitative morphology cut was derived 
for application to ECO. Traditional quantitative mor¬ 
phology discriminants, such as the concentration index 
C r = Rgn % /Rf,ci% defined using SPSS catalog photom - 
etry (e.g., Strateva et al. 2001 Shimasaku et al. 2001), 
yield unfortunately high error rates in the ECO sample 
(see Fig. [7p,j[^] We instead employ the pa metric recently 
developeaby K13, which combines the surface mass den¬ 
sity within i? 9 o% and a multiple of 1.7 times the differ¬ 
ence between the surface mass densities within f? 50 % and 
within the f? 50 % — Rgo% annulus. We optimize this met¬ 
ric for use as an early/late type discriminant in ECO 
by considering it as a function of M r and choosing two 
linear cuts in this parameter space that yields the min¬ 
imum misclassification rate in two separate M r regimes 
(see Fig. &)• Our optimized morphology discriminant is 


8 Note that we explicitly avoid the inclusion of color as a param¬ 
eter used for morphology discrimination due to the bias it would 
introduce against blue early and red late types. 


defined by: 


MA = 


-0.12 x M r + 7.1 M r < -19.5 
-0.49 x M r - 0.2 M r > -19.5. 


( 2 ) 


To reduce the misclassification rate resulting from im¬ 
plementing a quantitative morphology cut, we make fur¬ 
ther use of the Galaxy Zoo clean and debiased mor¬ 
phology classifications. Where such classifications ex¬ 
ist for ECO galaxies, we use the Galaxy Zoo early/late 
type classification rather than that inferred from the op¬ 
timized quantitative morphology cut described above. 
This substitution results in a naive apparent misclassi¬ 
fication rate of ~3.19% for late types and ~15.0% for 
early types in the morphology calibration sample. 

We can better estimate the misclassification errors that 
would result from applying this quantitative calibration 
to an independent sample with a bootstrap resampling 
approach applied to the calibration sample, usi ng th e 
“.632” error rate estimator as described by Efron (1983). 
The bootstrap sampling procedure is repeated for 1000 
iterations. For each iteration, we randomly sample N 
objects from the morphology calibration sample with re¬ 
placement, where N is equal to the total number of ob¬ 
jects in the calibration sample (the “.632” nomenclature 
refers to the fact that ~63.2% of the objects will end up 
in each bootstrap sample on average because objects can 
be selected twice). We then use each bootstrap sample 
to determine an optimum classification rule and evaluate 
the misclassification rate using this rule when applied 
to those members of the morphology calibration sample 
that were outside the bootstrap sample. As for the ECO 
sample, we include replacement of quantitatively inferred 
classifications with those from Galaxy Zoo to fairly as¬ 
sess the error rate during each iteration. We calculate 
our final error estimates as: 


E final 0.632 Eb avg + 0.368 E app , (3) 

where Eb aV g represents the average of the estimated er¬ 
ror rates over the 1000 resampling iterations and E app 
represents the naive apparent error rate estimate in the 
calibration sample. The final error estimates are similar 
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Fig. 8. — Illustration of the distributions of group membership 
that correspond to low and high estimated group halo mass regimes 
in our group finding analysis. 



Log( group N) 


Fig. 9.— Illustration of the recovery of real-space mock galaxy 
catalog group velocity dispersions after running group finding on 
a redshift-distorted version of the catalog and applying our finger- 
of-God collapse procedure on the resulting group catalog. Plot¬ 
ted are the distributions of real-space group velocity dispersion 
(red/orange scale) and corrected group velocity dispersion (green 
scale) versus real group N, where color coding of points represents 
the histogram density normalized to sum to one within each bin of 
group N (darkest colors imply the highest density). Red and green 
points are plotted offset by half the group N bin size for clarity. 
Note that we typically do not collapse groups to velocity dispersion 
levels as small as those seen in real-space groups. 

to the aforementioned naive apparent error rates, with 
estimated ~3.23% error rate for late types or ~15.8% 
for early types. We also calculate these error estimates 
in multiple individual bins of group halo mass and galaxy 
stellar mass as the error rate is not necessarily constant 
for different subclasses of galaxies in the sample. Our 
early-type misclassification rates tend to decrease with 
increasing galaxy stellar mass, from ~19% at our lowest 
stellar masses to ~11% at our highest stellar masses, and 
to increase with increasing group halo mass, from ~14% 
at our lowest halo masses to ~17% at our highest halo 
masses. 


3.5. Identification of UV Disks 

For the ECO galaxies covered by GALEX NUV imag¬ 
ing, we apply an automated procedure for quantitatively 
identifying Ultraviolet-Bright (UV-B) disks according to 


the definition of Moffett et al. (2012). For each galaxy, we 


apply the center position, position angle, and ellipticity 
of the SDSS optical ellipse fit to a set of fixed param¬ 
eter IRAF ellipse fits on the NUV imaging, which are 
allowed to proceed outwards radially until no significant 
UV flux is detected. We apply additional SExtractor- 
derived masking of the UV images beyond that deter¬ 
mined for the original SDSS optical photometry proce¬ 
dure, necessary due to the low resolution of the GALEX 
data and the occasional appearance of new contaminat- 
in g sources n ot p resent in the optical data. As described 
by Moffett et al. (20121, the quantitative UV-B disk clas¬ 
sification we employ requires satisfaction of an NUV— K 
color condition, which is included to ensure young stel¬ 
lar population ages and select stellar populations with 
a minimum ~10% young component by mass. We cal¬ 
culate NUV—A' colors between the optical g-band 50% 
light radius and the end of the NUV profile and require 
NUV—A' < 4.5 for classification as a UV-B dis!0 
Validating this fully aut omated identification a pproach 
against the methods of Moffett et al. (2012|, which 
employed a more detailed, galaxy-by-galaxy elliptical 
isophote fitting procedure, our new algorithm identifies 
up to ~25% more UV-B disks in an E/SO sample. Rea¬ 
sons for the difference between the more and less auto¬ 
mated approaches include occasional position angle mis¬ 
alignments between the galaxy regions where optical and 
UV emission dominate and imperfect automated mask¬ 
ing of UV sources, which is in some cases less aggre ssive 
than the masking employed by Moffett et al. (2012). 


3.6. Environment Metrics 

For the catalog samples under consideration, we com¬ 
pute two different metrics of galaxy environment. We 
primarily focus on environmental trends using group halo 
mass, but we also compare to results obtained with a 
smoothed galaxy density field in f|4] 

3.6.1. Group Finding 

We identify grou ps of galaxies usin g the friends-of- 
friends algorithm of Berlind et al. (2006). We infer group 
halo masses based on the observed total r-band lumi¬ 
nosity of galaxies in each group following an abundance 
match ing procedure, as described in |Blanton fe Berlind] 
(20071. The calculation of total luminosities for each 


group relies on high completeness to a fixed luminosity 
floor, so our group catalog is constructed for all ECO 
galaxies with M r < —17.33, even though we further an¬ 
alyze only a subsample of these galaxies satisfying our 
Mbary limit. Our group halo mass estimates assume a 
monotonic relation between the observed group luminos¬ 
ity and its halo mass. This relationship is determined 
by matching the space density of observed groups of a 
given luminosity to the theoretical space density of dark 
matter halos of a given mass, as derived from a standard 
concordance cosmology halo mass function. From consid¬ 
eration of the mock galaxy catalogs described below, we 
find that typical group halo mass errors are of order 0.12 
dex, although much larger errors can sometimes occur 


9 Note that in some cases galaxy outer disks are not well de¬ 
tected in our 2MASS images. In such cases, we calculate upper 
limit K-band magnitudes in the outer-disk regions and use these 
magnitudes for determining the NUV—A' color. 
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Fig. 10.— Illustration of the finger-of-God collapse procedure applied to the ECO catalog. Red points show the original inferred spatial 
distribution of galaxies, and smaller green points show the distribution after the finger-of-God collapse procedure has been applied. 


when groups are improperly fragmented or improperly 
linked together^] We label the brightest galaxy in each 
group the central and all other galaxies as satellites, and 
imperfect group identifications can affect these central 
vs. satellite designations, adding noise to comparisons 
of central vs. satellite trends. Note that in some cases 
the galaxy collections we refer to here as “groups” can 
in fact consist of single detected galaxies (that is, groups 
with N=l), which occurs most often for the lowest mass 
“groups” we probe (see Fig. [8]). 

The group-finding algorithm we employ automatically 
determines an appropriate “linking length” for grouping 
individual objects together, equal to b± and 6|| times the 
mean separation between galaxies in the input sample in 
the on-sky and line-of-sight directions, respectively (with 
D = cz/Hq distances and H 0 = 100 h km s -1 Mpc” 1 
us ed for this calculat i on). F ollowing the recent analysis 
of Duarte & Mamon (2014) we use values of b± = 0.07 


10 Note that the zero point of our group halo mass scale is de¬ 
pendent on comparisons with the mock catalog, as we calibrate our 
group halo masses to best approximate the true group halo masses 
in the mock catalog for the majority of halos in our sample, i.e., 
those at low mass. We subtract a simple constant correction of 
~0.15 dex from our initial abundance-matched mass estimates to 
correct for a tendency towards overestimation at low group halo 
mass. A significantly more complicated correction would be re¬ 
quired to yield perfect reproduction of the true group masses in 
all mass regimes, but using this simple correction we obtain me¬ 
dian estimated-mass-minus-true-mass offsets of only 0.02 dex for 
Mhalo < 10 12 Mg and -0.12 dex for Mhaio > 10 12 Mq, which are 
less than or consistent with the typical group halo mass errors 
overall. 


and fey = 1.1 for our group finding, as these values are 
found to be optimal for creating group catalogs that will 
be used for statistical study of the effects of galaxy en¬ 
vironment. For a sample region as small as RESOLVE- 
B, this mean spacing determination is far more sensi¬ 
tive to cosmic variance than for a large region as in 
the ECO sample. To compensate, when we apply the 
group-finding algorithm to RESOLVE-B we fix the link¬ 
ing lengths to those determined from a version of the 
ECO catalog limited to M r < —17, the magnitude limit 
of RESOLVE-B. We also use the luminosity to group halo 
mass conversion determined from this version of the ECO 
group catalog, since this conversion is similarly sensi¬ 
tive to cosmic variance. Through testing with an “ECO- 
analog” version of RESOLVE-B (matched to ECO’s shal¬ 
lower depth and lower completeness; see (3.8), we find 
that the difference in completeness between KESOLVE- 
B and ECO does not significantly bias the group halo 
masses we estimate, with scatter between the group halo 
masses inferred from the two versions of RESOLVE-B 
only reaching ^0.2 dex at Mhaio <1 1O 13 M 0 . The pri¬ 
mary effect of extra completeness on group identifica¬ 
tions is that additional N = 1 halos, typically containing 
faint galaxies near the magnitude limit, are found. Fig. 
[5] reveals that ECO is not complete for group halos with 
masses less than Mhaio ~ 10 11 Mq, so we refrain from 
including group halos below this mass limit in our anal¬ 
ysis. 


3.6.2. Smoothed Galaxy Density Field 
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We calculate a smoothed galaxy density field with 
an IDL p rocedu re based on the approach of Grogin| 
& Geller[ ([1999} . This procedure takes in individual 
galaxy redshift-space positions, assuming line of sight 
D = cz/H 0 for consistency with other methods applied 
here, and creates a continuous number density field from 
this spatial point distribution by smoothing each galaxy 
with a unit-normalized Gaussian kernel and summing the 
resulting galaxy space and velocity distributions. Since 
the samples considered in this work are volume- and 
not flux-limited, no luminosity function weighting fac¬ 
tors have^ been applie d in the density field calculation 
(cf. Grogin fc Geller| 1999;). 

In using this smoothed galaxy density field procedure, 
we use a smoothing kernel width of ~1.43 Mpc (1 Mpc/h 
with h = 0.7), a scale that is similar to a typical group 
halo virial radius at z ~0. In all subsequent density 
analysis, we limit ourselves to consideration of galaxies 
that lie >2 smoothing lengths from catalog edges and 
report density values normalized by the median of the 
density field values associated with these galaxies. 

In order for densities smoothed on this small scale to 
be physically meaningful, we also implement a proce¬ 
dure to statistically collapse “fingers of God” in ECO, 
as these cause measured line-of-sight velocities in large 
groups/clusters to imperfectly reflect physical distances. 
Similar methods for statistically correcting redshift-space 
distortions have bee n applied by other authors (e.g., 
Tegmark et al.||2004|. The algorithm we apply was de- 
veloped and calibrated on mock galaxy catalogs derived 
from an N-body simulation of a ACDM cosmological 
model, with Q m =0.25, Oa =0.75, fit, =0.04, h=0.7, 
n s =1.0, and ag =0.8. Initial conditions were s et using 


seco nd order Lagrangian Perturbation Theory (Scocci- 


marro 


1998) at a starting redshift of z =99 and the par- 


ticfe distribution was evolved using the code GADGET- 
2 (Springel 2005). The simulation contains 1050 3 dark 
matter particles in a box of size 180 Mpc/h, sufficient to 
encompass multiple ECO catalog volumes. The result¬ 
ing particle mass is 3.5x 10 8 Mq /h and the gravitational 
force softening is 7 kpc/h. At this resolution, the low¬ 
est mass halos in the ECO catalog sample typically con¬ 
tain 160 particles. We identify halos in the dark matter 
particle distribution using a friends-of-friends algorithm 
with a linking length equal to 0.2 times the mean inter¬ 
particle separation. We then populate these halos with 


gal axies using a halo occupation distribution (Berlind & 
Wein berg|2002 ) designed to produce a galaxy population 
with - smnlar space density and clustering properties as 
ECO galaxies. Central galaxies are given positions and 
velocities of halo centers of mass, and satellite galaxies 
are given positions and velocities of randomly selected 
dark matter halo particles. This procedure produces a 
“real-space” mock galaxy catalog. The real-space cata¬ 
log is subsequently distorted into a ’’redshift-space” mock 
galaxy catalog by assuming the line of sight direction to 
extend radially outward from the center of the box and 
incorporating galaxy peculiar velocity components along 
the radial direction into galaxy redshifts. 

Comparisons between real-space and redshift-space 
versions of the mock galaxy catalogs allow us to deter¬ 
mine the distribution of 3D, real space radius offsets from 
the group center positions as a function of the group 
virial radii ( R V i r ), with R V i r estimated using the group 


masses and the defining group overdensity. The over¬ 
all real-space group-centric R/R v t r density profile deter¬ 
mined from the mock catalogs is approximated with a 
gamma probability density function fit and used as an 
input when we apply the collapse algorithm on observed 
catalogs. The fit relation has the following form: 


f(x,a,b) = 


x a 1 x e bx x b a 


T(a) 


(4) 


where a=1.6 and b=5. 

We begin the finger-of-God collapse process on the ob¬ 
served galaxy catalog s with group identifications deter¬ 
mined as described in §3.6. 1| For each galaxy in an N > 1 
group, we first assign a random real-space displacement 
from the center of the group according to our fit distri¬ 
bution of mock-catalog-determined R/R V i r values, tak¬ 
ing into account the group’s estimated R V i r . Each as¬ 
signed 3D radius is then checked for consistency with 
the galaxy’s observed projected distance from the group 
center, that is, a galaxy’s assigned 3D distance from the 
group center is required to be greater than its observed 
projected distance from the center. If this is not the 
case on first assignment, a different random 3D radius 
value is chosen until the condition is met. This proce¬ 
dure does not significantly bias the assigned R/R vir dis¬ 
tribution. Next, the appropriate redshift direction dis¬ 
placement from the group center is determined for each 
galaxy, using the observed spatial coordinates and the as¬ 
signed 3D radius. The sign of the redshift displacement 
for each galaxy relative to the group center is determined 
randomly. 

Validation of the final collapse procedure on the mock 
catalogs yields reasonable agreement between the cor¬ 
rected and original simulated, real-space group velocity 
dispersions (see Fig. [9]), although we tend not to reach 
corrected group velocity dispersions as small as the orig¬ 
inal dispersions seen in real-space mock catalog groups). 
This bias appears to be due largely to imperfect recovery 
of group halo R V i r by our group finding method; over¬ 
estimation of R v ir has the effect that our collapse al¬ 
gorithm places galaxies at larger 3D group-centric radii 
(and larger velocities) than observed in the real-space 
groups. The effect of the collapse process on the distri¬ 
bution of galaxies in the ECO catalog is shown in Fig. 
10 The application of this collapse process yields greater 
confidence in our density field results, which generally 
agree with group-finding results that are less sensitive to 
such finger-of-God effects. 

3.7. Atomic Gas Mass Estimates 


We cross-m atch to published ALFALFA a.40 (Haynes] 
et al. 2011) HI sources to find possible ECO galaxy 
matches and check for possible confusion usin g a large 
combined redshift catalog containing ECO (see 1 2.1). We 
flag a source as confused if more than one galaxy is found 
within the 3' matching radius (1.5 times the effective AL¬ 
FALFA beam radius) and if the HI profiles would over¬ 
lap in velocity space assuming minimum 50 km/s red¬ 
shift uncertainties and typical line widths of ~100 km/s 
(the weighted aver age of the ALFALFA-d erived veloc¬ 
ity width function; Papastergis et al.||2011). We use the 
reported line flux densities to convert to Mhi via the 
standard formula ( Mhi = 2.36 x 10 5 x D 2 x Fhi Mq; 
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Fig. 11. — Multiplicative completeness correction factors derived for each individual ECO mass-limited sample galaxy (points) as a 
function of r-band absolute magnitude, model g — i color, and r-band surface brightness ^ 23 . 75 r- Galaxies indicated by red points require 
further multiplicative correction factors on top of the factors shown, due to the loss of galaxies with extreme line-of-sight velocities outside 
the ECO definition. The majority of galaxies indicated by red points lie in the Coma Cluster, where the additional multiplicative correction 
factor reaches its maximum value of ~ 1.42. 


Haynes & Giovanelli 1984). We then multiply our masses 
by 1.4 to correct tor the presence of He. In the absence of 
HI detections, we estimate upper limits according to the 
procedure described by K13, which uses the typical a.40 
rms noise as a function of declination integrated over 
a velocity interval estimated for each source according 
to the average relation between internal velocity and r- 
band magnitude (logV = —0.29 — 0.123 M r , calibrated 
by K13). 

Where ECO galaxies lack ALFALFA a.40 detections 
or have confused HI measurements, we use the “p ho- 


tometric gas fraction” technique (Kannappan 2004) to 
assign HI mass estimates based on the observed tight 
correspondence between HI gas to stellar mass ratio and 
u — J color. We employ the photometric gas fraction cal¬ 
ibration and procedure described by K13, in which gas 
masses are assigned to galaxies with (it — j) modei < 3.7 
according to the relation: 


log (M ff j/M.) = 2.7 - 0.98 (it - J) 


model 


(5) 


with random 0.34 dex scatter motivated by the observed 
relation. If the gas mass estimated with this relation for a 
given galaxy would exceed its calculated upper limit, the 
upper limit is adopted instead. For galaxies with colors 
redder than (it — j^rnodei _ 3 ^ that ^ w j lere the linear 
color-gas-to-stellar mass ratio relation breaks down for 
quenched galaxies, the procedure assigns random values 
in the logarithmic range Mhi/M * = 0.001 — 0.5, again 
constrained to lie below the estimated upper limit for 
each galaxy. 

We make one significant modification to the K13 pho¬ 
tometric gas fraction procedure, which is motivated by 
the tendency for the bluest and most gas-rich galaxies 
to lie above the aforementioned gas-to-stellar mass ra¬ 
tio versus color relation (see K13 Fig. 8 a). To better 
reflect the typical gas-to-stellar mass ratios observed for 
galaxies with (u — J) modei < 3 ? we multiply the gas mass 
estimates that would result from the K13 estimator by an 
additional factor of 1.5 in this regime, with resulting es¬ 
timates still constrained to lie below upper limits where 
available. While our gas mass estimates are improved 


by the inclusion of this factor, we note that its inclusion 
does not qualitatively affect the results we report, which 
remain similar even if the extra multiplicative factor is 
omitted from estimated gas masses. In particular, we 
largely focus our analysis of gas content on the incidence 
of gas-dominated (Mhi/M* > 1) galaxies in ECO, and 
so for much of our analysis, the exact gas mass estimated 
for a given galaxy is less important than whether or not 
it falls into a broad gas-to-stellar mass ratio category. 

3.8. Completeness Corrections 

Since the RESOLVE-B sample has much higher com¬ 
pleteness than the ECO sample (see j ]2.2[ ), we use the 
former to calibrate and correct for the effects of redshift 
incompleteness in the latter. We do so in two steps: 
first we compare the full RESOLVE-B galaxy number 
counts with a version of RESOLVE-B that contains only 
galaxies with DR7 redshifts (that is excluding the ex¬ 
tra SDSS Stripe 82 redshift coverage beyond the original 
redshift survey), and next we carry out a similar com¬ 
parison between ECO catalog galaxy number counts and 
those of a version of ECO that contains only galaxies 
with DR7 redshifts. By dividing the RESOLVE-B DR7- 
to-full catalog correction factor by the ECO DR7-to-full 
catalog correction factor, we obtain a measure of the fac¬ 
tor needed to correct the ECO catalog membership to 
the same completeness as the RESOLVE-B sample. For 
both steps of this process, we divide both samples into 
cells in two-dimensional parameter space grids (e.g., M r 
vs. color), using grids determined with a simple adaptive 
approach, which begins with subdividing each axis into 
four large cells. If more than a set minimum number of 
galaxies are present in a given cell (minimum is five for 
RESOLVE-B and 100 for ECO), the cell is subdivided in 
half iteratively until no further subdivisions are allowed 
by the minimum number of galaxies per bin condition. 
We interpolate each irregularly gridded dataset into a 
smoothly varying number density field in order to derive 
the final completeness correction fields. We have exam¬ 
ined multiple different 2D parameter space options for 
deriving these corrections and find that in all cases the 
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Fig. 12.— Illustration of ECO catalog galaxy distributions in the group halo mass and galaxy baryonic mass space, where early types of 
red/blue and central/satellite classes are shown on panel a, and late types of red/blue and central/satellite classes are shown on panel b. 
No catalog membership completeness corrections have been applied in this figure. 


resulting corrections yield overall similar results. In or- 
der to reduce noise in the final correction caused by low- 
occupation bins, we choose to average the corrections de¬ 
rived for individual galaxies using the M r vs. (g — j) modei 
and M r vs. /-i 23 . 75 r parameter spaces (where /i 23 . 75 r rep¬ 
resents an average r-band surface brightness evaluated 
within an outer galaxy ellipse corresponding to a 23.75 
magarcsec -2 surface brightness level, the lowest-surface- 
brightness-level ellipse we find reliably determined). 

Our completeness correction method results in multi¬ 
plicative correction factors that vary as a function of ab¬ 
solute magnitude, color, and ZJ. 23 . 75 r as illustrated in Fig. 


11 The median correction factor applied to a galaxy in 


our sample is ~ 1 %, but correction factors can reach up 
to ~ 68 %. Regions where correction factors are enhanced 
in this figure correspond to grid cells where genuinely 
large differences in the number counts of the comparison 
samples occur. Applying the derived completeness cor¬ 
rection factors to ECO results in galaxy luminosity and 
color distributions that are in typically good agreement 
with RESOLVE-B sample properties. In addition to the 
general redshift catalog incompleteness, we note that the 
presence of cluster fingers of God in the ECO sample can 
cause further incompleteness through the loss of mem¬ 
bers with extreme line-of-sight velocities, putting them 
outside our sample limits even though we include a large 
buffer region in our analysis to mitigate such losses. The 
presence of the Coma Cluster, in particular, leads to in¬ 
completeness for the highest group halo mass in the ECO 
sample (see upper leftmost panel in Fig. [2j|, which can¬ 


not be appropriately corrected through consideration of 
the RESOLVE-B comparison sample as this incomplete¬ 
ness is caused by ECO boundaries and not incomplete 
redshift data. 

To quantify correction factors for this boundary effect, 
we construct a comparison catalog including the same 
data sources as ECO but extending beyond the ECO 
redshift limits in both directions (1500 km/s to 12000 
km/s), encompassing all apparent fingers of God extend¬ 
ing outside of the ECO region. As reprocessed photom¬ 
etry is not available for this entire comparison catalog, 
we rely on SDSS catalog magnitudes and perform group 
finding on the comparison catalog using catalog r-band 
magnitudes. We fix the linking lengths for this group 
finding to the values appropriate for the normal ECO 
region to avoid adverse effects from the fall-off in com¬ 
pleteness beyond the ECO boundaries. For any galaxy in 
a group with N>10 that has been affected by proximity 
to the ECO volume edges, we calculate and apply an ad¬ 
ditional completeness correction factor based on the ratio 
of the number of galaxies in the full comparison catalog 
group to the number of galaxies within the ECO cat¬ 
alog boundaries, derived considering only galaxies with 
catalog magnitudes M r < —18.4 (corresponding to the 
nominal SDSS r-band completeness limit at cz = 12000 
km/s). The boundary completeness correction factors 
are only necessary for two groups, with a maximum cor¬ 
rection value of ~ 1.42 for Coma and a correction value 
of ~ 1.14 for the NGC4065 group (see Fig. [ 2 ]). We note 
that these correction values are of necessity approximate 
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Fig. 13. — Illustration of the traditional morphology-environment relation in the ECO and RESOLVE-B samples. Early and late type 
frequencies as a function of group halo mass in (a) ECO and (b) RESOLVE-B, with early and late type frequencies crossing over at 
Mhalo ~ 10 13 5 Mq in the ECO sample. Greyscale dotted lines indicate the frequencies for central galaxies alone, which become poorly 
determined at high group halo masses due to the small number of high halo mass groups (and thus centrals) present in the sample. In panel 
a, all frequencies are plotte d at their “expected” value given the calibrated uncertainties in our semi-quantitative morphology classification 
method, described in §3.4| The error bars in this panel are plotted as a combination of the calculated misclassification errors in each group 
halo mass regime and the (binomial) counting statistics in each bin, while those in panel b reflect only the relevant counting statistics 
(assuming zero misclassification). Note that the early/late type frequencies for the two samples do not strictly agree in all bins, particularly 
near Mh a ; 0 ~ 10 12 Mq. Applying the quantitative classification approach to RESOLVE-B results in frequencies that more closely agree with 
the frequencies in ECO, although slight residual differences remain, perhaps resulting from group-to-group variations in typical properties. 
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Fig. 14.- Illustration of morphology-environment relations in 
the ECO sample for separate high-mass ( Mf, ary > 10 1 0 A4q ) and 
low-mass (M^ ary < 10 1o Mq) galaxy subsamples, indicated by 
thick and thin lines, respectively. Since few high mass galaxies 
inhabit group halos below ~ 10 12 Mq, we do not plot the frequen¬ 
cies for high mass galaxies below this point. While the trends are 
qualitatively similar over most of the group halo mass range, the 
relations have different amplitudes for the low and high galaxy 
mass samples. 


as they do not account for any possible population differ¬ 
ences between the parts of a cluster interior and exterior 
to our sample volume. As they apply to only two group 
halos with M halo > lO 14 Af 0 in our sample, a regime 
in which it is already difficult for us to draw conclusions 
due to the small number of rich clusters sampled in a vol¬ 
ume of ECO’s size, we caution against over-interpreting 
results solely from this cluster regime. We will largely 
focus on more common moderate to low richness envi¬ 
ronments in this work. 

As the loss of galaxies beyond ECO boundaries affects 
sample membership, this loss can also affect our derived 
group halo mass estimates. To correct for this effect, we 
use the same logic described above and apply a correc¬ 


tion factor to each affected group halo mass based on the 
ratio of the total group luminosity derived in the com¬ 
parison catalog to the total group luminosity derived in 
the ECO catalog. The applied group halo mass correc¬ 
tion for Coma is an increase of ~ 0.1 dex but is negligible 
for the NGC4065 group. 


4. RESULTS 

In this section, we explore trends in galaxy properties 
related to galaxy mass and environment with an eye to¬ 
wards their connections to galaxy disk (re)growth. In all 
plots of our results we include only objects that meet our 
final group membership and Mt, ary > 10 9 ' 3 M 0 selection 
as detailed in fj2l We use group halo mass as our pri¬ 
mary indicato r of en vironmental richness. As mentioned 
previously in i 3.6.1 we also consider isolated galaxies as 
members of groups with N=l. Although this description 
may at first seem counterintuitive, it does provide a way 
to quantify group environment as a continuous variable 
through group halo mass, which is useful as a basis for 
comparison with theory. Since the r-band luminosity of a 
galaxy provides an excellent proxy for its total baryonic 
rather than stellar mass (as discussed further in Kan- 


nappan et al. 2013), it is the case that where group IN=1 
our estimated group halo mass is directly related to the 
lone galaxy’s baryonic mass. We can nonetheless identify 
group halo environment rather than galaxy mass as the 
primary driver of a trend if it is sensitive to group central 
versus satellite status at fixed galaxy mass. 

The distribution of various classes of ECO catalog 
galaxies in group halo mass versus galaxy baryonic (stel¬ 
lar plus atomic gas) mass parameter space can be seen 
in Fig. [12] Results derived from the density field are 
qualitatively similar to results derived from group halo 
masses, except as discussed in this section. 


4.1. Traditional Morvholoqy-Environment Relation - 

P(M\E) 

























































































































14 


Moffett et al. 


o 

c 

CD 

3 

cr 

CD 



9.5 10.0 10.5 

Log( M bary [M 


11.0 


sun. 


]) 


9.5 10.0 10.5 

Log( M bary [M 


11.0 


sun, 


]) 


Fig. 15.— Illustration of morphology-mass relations in the ECO sample for multiple color/morphology subclasses of ECO galaxies, where 
the left panel includes centrals only and the right panel includes satellites only. 


As seen in Fig. for the traditional formulation 
of the morphology-environment relation (frequency of 
a particular morphology as a function of environmen¬ 
tal richness), the ECO sample displays the expected in¬ 
crease in early-type and decrease in late-type frequen¬ 
cies as a function of increasing environmental richness, 
here represented by increasing group halo mass. Near 
group halo mass ~ 10 13 ' 5 Mq, we observe a crossover 
point where early and late type frequencies become ap¬ 
proximately equal. The frequencies in the more complete 
RESOLVE-B sample, where morphological classifications 
are entirely based on by-eye judgments, do not strictly 
agree with the frequencies in E CO, particularly around 
Mhaio ~ 1O 12 M 0 (see Fig. 13 d). If we apply quantita¬ 
tive classification methods to RESOLVE-B, we find fre¬ 
quencies that more closely but not completely agree with 
the frequencies in ECO. Variations in morphological mix¬ 
tures between groups at fixed halo mass could also plau¬ 
sibly contribute to differing average morphological mixes 
in these two samples. If we divide central and satellite 
galaxies, we find that central and satellite morphology- 
environment trends are similar to each other in both sam¬ 
ples, where the satellite galaxy tre nd ( not shown) closely 
follows the combined trend in Fig. 13i/b. If we examine 


the traditional morphology-environment relation for low 
and high baryonic mass (divided at Mbary = 1O 1O M 0 ) 
galaxies separately, we qualitatively recover the expected 
early-/late-type frequency trends as a function of group 
halo mass, with the high baryonic mass relations shal¬ 
lower than the low baryonic mass relations and offset to¬ 
wards higher early-type frequency. The early/late-type 
frequencies cross over at group halo masses > 10 13 Mr. 


01 


and both trends steepen in this regime (see Fig. 14). 

Dissecting the morphology-environment and 
morphology-mass relations furt her by including galaxy 
color, we can see from Fig. [15] that the offset between 
low and high baryonic mass morphology-environment 
relations must be partially driven by the significant 
changes in blue and red galaxy frequencies as a function 
of galaxy mass, specifically as blue late types are 
significantly more common at low galaxy mass and red 
early and late types are more common at high galaxy 


mass. In Figs. |T5] and |T6j we also see that centrals and 
satellites within each galaxy color/morphology class 
follow qualitatively similar trends as a function of galaxy 
mass and environment. However several nota ble offsets 
in frequency occur, for example in Fig{T5j blue late 
types are less common among satellites than centrals 
at low baryonic mass, whereas red early types are more 
common among satellites than centrals at low bar yonic 
mass. Interestingly, in the top panels of Fig. [l6j blue 
late types are more common among satellites than 
centrals at low group halo mass, illustrating that galaxy 
demographic trends in the low galaxy mass and low 
group halo mass regimes are not necessarily equivalent. 
There is also a hint of a similar reversal for red early 
types (i.e., red early types are marginally more common 
among centrals than satellites at low group halo mass). 
As can be seen in the bottom two panel sets in Fig. 
16 the galaxy color/morphology class frequency trends 
as a function of group halo mass, where they can be 
quantified, are overally similar between low and high 
galaxy mass subsamples. 

4.1.1. Focus on Blue-sequence Early Type Galaxies 


As illustrated in Fig. [TTp ,, the majority of late-type 
galaxies occupy Mhaio < 1O 12 M 0 environments, which 
is not surprising given the form of the traditional 
morphology-environment relation. More interestingly 
the same figure shows that ~50% of early-type galax¬ 
ies occupy Mhaio < 10 12 Mq environments. As seen in 
Fig. [IT]), many early types in the lowest density environ¬ 
ments are blue-sequence early types. These blue early 
types become most common below M ha i a ~ 1O 11 ' 5 M 0 . 
Blue early types with baryonic masses large enough to 
meet our mass limit are primarily central galaxies in this 
low group halo mass regime and primarily satellite galax¬ 
ies in richer environments (see Fig. 12). However, since 
satellite blue early types with masses below our survey 
limit could also populate the low mass halos, we cannot 
yet quantify the balance between blue early type cen¬ 
trals and satellites in these environments. Fig. |T8| shows 
that blue early-type galaxies appear to have similar envi¬ 
ronment distributions to blue late-type galaxies overall, 
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Fig. 16. — Illustration of morphology-environment relations in the ECO sample for multiple subclasses of ECO galaxies. All plots in the left 
column include centrals only, while plots in the right column include satellites only. The bottom two panels show morphology-environment 
relations in the ECO sample for separate low-mass (Mfc ar?/ < 1O 1O M0) and high-mass (Mfc ar2/ > 1O 1O M0) galaxy subsamples. 
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but blue early types display a slight preference towards 
lower group halo mass environments. As calculated with 
the Kolmogorov-Smirnov (K-S) test, the group halo mass 
distributions of blue early and blue late types in Fig. [18] 
have ~ 1% probability of being drawn from the same 
parent distribution (P same ~ 0.01). 

We find that blue-sequence early-type galaxies are 
not only most common in low group halo mass envi¬ 
ronments but at low galaxy masses as well, both of 
which are regimes where extreme gas richness is typi¬ 
cal. As was previously found by KGB, we observe that 
the blue early-type galaxies in ECO become more com¬ 
mon with decreasing stellar or baryonic mass (Figs. 19 
and 16 1 . These galaxies only emerge in ECO with mea¬ 
surable frequency around the galaxy bimodality mass 
(Mb ~ 10 1o 5 Mq) and increase in frequency significantly 
below M* ~ 10 10 Mq. This behavior is similar to that 
observed in KGB, where the blue early-type population 
increases sharply below M* ~ 10 9 ' 7 M 0 , the “gas-richness 
threshold” mass. Our observed frequency transition is 
somewhat less sharp than that observed by KGB, but it 
is plausible that any sharp transitions could be washed 
out by the error rates inherent in our semi-quantitative 
classification method. Our absolute frequencies are sim¬ 
ilar to those reported by KGB, although not in strict 
agreement within our error bars. If similarly large error 
bars on the KGB frequencie s are assumed, then th e two 
trends would agree overall. Thomas et al. (2010) have 
also observed an increasing frequency ot “rejuvenated” 
early-type galaxies with decreasing galaxy mass (and en¬ 
vironmental density), reaching a maximum of ~45% of 
the early-types, which is similar to our observed maxi¬ 
mum frequency. This observed low-mass preference im¬ 
plies that mass-dependent mechanisms are closely tied 
to the rise of the blue-sequence early-type population. 
We note that the typically low stellar mass of blue early- 
type galaxies appears linked to the typically low group 
halo mass environments they inhabit. At M* ~ 1 O 19 M 0 , 
the typical baryonic mass for a galaxy is ~ 10 10 ' 1 to 


traditional morphology-environment relation is approx¬ 
imately flat until group halo masses above ~ 10 12 5 M 0 


10 w ' 2 Mq, which from Fig. 12 corresponds to a typical 


group halo mass of ~ 1 O 115 M 0 for centrals. 

4.2. Alternative Morpholoqy-Environment Relation - 

P(E\M) 

An alternative formulation of the morphology- 
environment relation, expressing the probability for a 
galaxy with a given morphology to inhabit a particular 
environment, is illuminating for understanding the “typ¬ 
ical” environments of various galaxy classes as a function 


20 } As seen in this 
ue early types and 


of galaxy mass as illustrated in Fig. 
figure the typical environments of bl 
blue late types are similar at constant baryonic mass, for 
both centrals and satellites. This figure also illustrates 
the general lack of a P(E|M) morphology-environment 
relationship for central galaxies at fixed mass, which also 
holds if we consider all early types and all late types to¬ 
gether and is essentially a manifestation of the known 
central galaxy mass to group halo mass relation (e.g., 
Yang et al. 2009 Most or et al.||2010 ). 

The traditional F(M | E) morphology-environment rela¬ 
tion obscures this fact by mixing together morphological 
classes with differing baryonic mass distributions within 
each environment regime. However, for low baryonic 
mass (Mbary < 1 O 1 O M 0 ) galaxies taken alone even the 


(see Fig. 14), about 1 dex into the regime where these low 
mass galaxies have become satellites. Accordingly, in the 
P(E|M) formulation, a morphology-environment relation 
reemerges for low baryonic mass satellites: in general 
early-type satellites typically occupy higher group halo 
mass environments than late-type satellites at fixed bary¬ 
onic mass, although specifically this relation is driven by 
red early types. For satellite galaxies, the strongest dif¬ 
ference in typical group halo mass is between red galaxies 
and blue galaxies, where both red early and late types 
occupy higher group halo mass environments than either 
blue early or late types. 

Even though low mass blue early and blue late types 
occupy environments that are typically similar (and are 
classified as centrals with the same ~80% frequency for 
Afhary < 1 O 1 O M 0 ), their full environment probability dis¬ 
tributions at a given mass may not necessarily be the 
same. In this case, we find K-S test P sa me ~ 0.22 for 
the group halo mass distributions of Mbary < 1O 1 O M 0 
blue early and blue late types, where the blue late type 
environment distribution appears to be broader, suggest¬ 
ing greater environmental diversity at fixed mass (al¬ 
though we note that the blue late type distribution is 
better sampled). Considering environmental density val¬ 
ues (smoothed on ^1.43 Mpc scales), we find that the 
typical densities around Mb ary < 1 O 1 O M 0 blue early- 

S e centrals are lower than those of late types (see Fig. 

). If number density is more sensitive to major merg¬ 
ers than group halo mass, then this tendency towards 
somewhat lower environmental densities could be consis¬ 
tent with blue early-type centrals existing as post-merger 
objects with their number of neighbors reduced by merg¬ 
ing. As with group halo masses, blue early-type and blue 
late-type satellites have similar typical densities in the 
Mbary < 1 O 1 O M 0 regime, but the overall blue early and 
blue late type density distributions in this mass regime 
are more plausibly different, with a P sa me ~ 0.06 as 
quantified with the K-S test. 

4.3. Extreme Gas Richness and Environment 

Considering ECO+A galaxies of all types together, 
Fig. [22] shows that the fraction of gas-dominated galax¬ 
ies (i.e., those with Mhi/M * > 1) is a strong function of 
environment in general. The shape of the trend differs 
as a function of group halo mass compared to smoothed 
density field values: the halo mass relation displays a rel¬ 
atively smooth rise in gas-dominated galaxy frequency 
with a steep rise below ~ 1O 114 M 0 , whereas the density 
relation shows an overall smooth rise in gas-dominated 
galaxies toward lower densities (compare two trends in 
Fig. 22 1 . The sharp increase in the frequency of extreme 


gas richness in our sample below Mhalo ~ 1O U ' 4 M 0 
marks this as a regime where fractionally large gas reser¬ 
voirs become a common feature of galaxies. Furthermore, 
we see that central galaxies are primarily responsible for 
the sharp increase in gas-donrinated galaxy frequency 
with group halo mass, wit hin the baryonic mass limits 
of our sample (see Fig. 23). 

If we consider the typical Mhi/M. * values as a func¬ 
tion of environment for individual galaxy types in the full 
ECO sample, we find that in low group halo mass envi- 
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Fig. 17.— Group halo mass distribution for early and late type galaxies in the ECO sample, with corresponding grey dashed lines indicating 
the distribution for group central galaxies of each type. Panel a shows all early and all late types together, and panel b breaks down the 
early types further by red and blue sequence membership. Approximately 50% of ECO early type galaxies occupy M^aio < 10 12 Mq 
environments, approximately 20% of which are blue early types. 
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Fig. 18.— Group halo mass distribution for blue early-type and 
blue late-type galaxies in the ECO sample, with corresponding grey 
dashed lines indicating the distribution for group central galaxies 
of each type. 


ronments blue early types display Mhi/M * values that 
are comparabl e or somewhat lower than those of blue late 
types (see Fig. 241. We also observe that gas-dominated 
blue early types, with HI gas-to-stellar mass ratios > 1, 
are found almost enti rely in low group halo mass envi¬ 
ronments (see Fig. 25), becoming most common below a 
group halo mass of ~ 10 11,5 Mg, where they are primar¬ 
ily group centrals (although again, any satellites would 
likely fall bel ow our mass limit). 

From Fig. [24j we find that satellite galaxies display 
typically higher gas-to-stellar mass ratios than centrals 
at fixed morphology, color type, and group halo mass. 
It is somewhat surprising that higher satellite gas con¬ 
tent appears to persist in rich environments, except for 
red early types where the central and satellite values are 
indistinguishable. Satellites are typically expected to be 
quenched, but we remind the reader that many of our HI 
gas masses are estimated using optical colors, and this 
photometric gas fraction techniq ue m ay be unreliable in 
dense environments (Cortese et al. 2011). Smoothed den¬ 
sity field trends are similar to but we aker than the group 
halo mass trends illustrated in Fig. 24 As illustrated 
in Fig. [261 these environmental effects on galaxy gas 


Fig. 19.— Variation in blue early type galaxy frequency as a 
function of stellar mass in the ECO sample. The solid line indi¬ 
cates the frequency of blue early types as a fraction of early types 
only, while the dotted line indicates their frequency as a fraction of 
all galaxy types. Frequencies are plotted at their expected values 
given the calibrated uncertainties in our s emi- quantitative mor¬ 
phology classification method, described in §3.4| Error bars shown 
are a combination of the estimated misclassmcation errors and the 
(binomial) counting statistics in each bin. 

fractions are not solely due to the connection between 
galaxy mass and group halo mass; while high and low 
gas fraction centrals display similar group halo masses at 
fixed baryonic mass, high gas fraction satellites inhabit 
significantly lower mass groups than lower gas fraction 
satellites. Thus for satellites, group halo mass is a more 
important driver of gas richness than galaxy mass. 


4.4. UV Disk Growth and Environment 
In the ECO+G subsample, we find that UV-B disks 
are relatively common (Fig- 27), occurring with 34^1% 
overall frequency. Fig. [27] shows that the frequency of 
UV-B disks in early types is greatest at low group halo 
masses. UV-B disk frequencies among early type galax¬ 
ies only reach measurable values in our sample below 
Mhaio ~ 10 12 M 0 , approximately the same group halo 
mass scale below which blue-sequence early types and 
gas-dominated galaxies emerge. Correspondingly, we 
find that UV-B disks are most common in early types 
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Fig. 20. — Characteristic distribution of group halo mass as a function of baryonic mass for different galaxy types. Solid lines indicate 
the running median group halo masses for centrals of each galaxy type, while dashed lines indicate the medians for satellites of each galaxy 
type. Error bars for each median point are estimated from the dispersion in properties in each bin. The background greyscale levels indicate 
the probability of inhabiting a particular group halo mass at a given baryonic mass, as the histogram densities used to set the greyscale 
have been normalized to one in each baryonic mass bin (darkest points imply the highest probabilities). 


at low galaxy mass. The ~ 30% frequency of UV-B 
disks we find among early types at low baryonic masses 
(see Fig. 28) is consistent with the reported 42+®% UV- 
B disk frequencies of Moffett et al. (2012) for early- 
type galaxies in a low-to-mtermedlate stellar mass sam¬ 
ple. Within the mass ranges we probe, UV-B disk hosts 
below Mhaio ~ 1O 12 M 0 are typically central galaxies, 
but additional satellite UV-B hosts with lower masses 
could also exist in these environments. Early type galax¬ 
ies hosting UV-B disks show a strong tendency to host 
larger HI gas reservoirs than galaxies without UV-B disks 
(K-S test Psame ~ 9 x lCT+k see Fig . 27^ 


ment with results from Moffett et al. (2012 f 


m agree- 
br low-to- 


intermediate mass early type galaxies! Overlapping in 
their typically gas-rich, low group and galaxy mass na- 


11 Note that this result also holds with similarly high significance 
if we restrict only to galaxies in the ECO+A subsample. 


ture, blue-sequence early-type galaxies are often UV-B 
disk hosts, with 74+®% of blue-sequence E/SOs in the 
full ECO+G subsample hosting UV-B disks, or 77+g% 
in the low baryonic mass regime below ~ 10 1o Mq. 

5. DISCUSSION 

In this section, we compare our results on relationships 
between galaxy properties and environments to previous 
results in the literature and to galaxy evolution scenarios. 

5.1. Morphology-Environment Relations 
5.1.1. Comparisons to Previous Results 

In general, we find that in the traditional P(M|E) 
morphology-environment relation formulation, our mea¬ 
sured early and late type frequencies behave in a man¬ 
ner similar to that observed in previous studies of this 
relation, for example, with late-type galaxy frequencies 
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Fig. 21.— Characteristic distribution of environment al d ensity (<■>-'1.43 Mpc smoothing kernel) as a function of baryonic mass for different 
galaxy types, with symbols and lines analogous to Fig. |20| 


decreasing from ~80% in the least rich environments to 



Dressleri|1980 

Postman & Geller 

1984 

Whitmore 

et al. 

1993). 

Examining early-/late-type frequencies as 


et al.| (|2009| and [Hoyle et al.| (p 
constant frequencies in the high 


2012) find approximately 


group halo mass regime 
(> 10 13 Mq), which we do not observe in the ECO sam¬ 
ple, even if we select only high baryonic mass galaxies to 
best compare with these studies (although the error bars 
are large as discussed further below). In addition to their 
default es timates for group mas s based on virial radius 


measures, Bamford et al. (2009) specifically test the use 
of summed luminosities as a proxy for group mass, us¬ 
ing a definition close to that we employ, and find a very 
slightly more pronounced trend in the early-type frac¬ 
tion with this proxy, alt hough still weaker than our full 
sample trend. Similarly, Poggianti et al. (2009) find no 
significant frequency trend with cluster velocity disper¬ 
sion but do find a trend with another proxy for group 


halo mass, X-ray luminosity. The typical galaxy masses 
considered by these authors are higher than those consid¬ 
ered in ECO, and when restricting to M bary > 10 lo M Q 
galaxies alone, our early-/late-type frequency trends at 
the highest group halo masses are relatively weak given 
the large error bars, e xcept in the high est group halo 
mass bin (see Fig. 14). Calvi et al. (2012) also find vari¬ 
ations with galaxy mass, with their intermediate stellar 
mass galaxies showing similar morphological mixes in all 
environments except in the most massive clusters, which 
is compatible with the behavior of our M bary > 10 1o Mq 
subsample. Bamford et al. (2009) likewise report that the 
form of the morphology-environment relation is strongly 
dependent on the stellar masses of the galaxies consid¬ 
ered, with shifts in the overall frequency levels between 
subsamples. We find trends similar to the Bamford et 
al. results, where morphology-environment relations are 
similar in shape but offset between low and high mass 
subsamples in the sense that late-type frequencies ar e 
typically higher among low-mass galaxies (see Fig. 14). 
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Fig. 22. — The environment-dependent frequency of HI gas-to- 
stellar mass ratios greater than one for ECO+A galaxies. The 
vast majority of ALFALFA a.40 derived HI limit values for ECO 
galaxies are sufficiently strong to place these galaxies in either the 
Mhi/M* > 1 or Mhi /M* < 1 category, but for any limit val¬ 
ues found in the Mhi/M * > 1 category, we replace th ese values 
with photometric gas fraction estimates as described in ^3.7| If we 
leave the limit values in place and use no photometric gas frac¬ 
tion estimates, we still find trend lines that are within the error 
bars of the trends shown here. Confused sources in the ECO+A 
sample have been omitted. The thick dark blue line illustrates the 
sharp uptick of the gas-dominated galaxy frequency at low group 
halo mass (lower x axis), and the thin light blue line illustrates the 
smoother increase in gas-dominated galaxy frequency as a function 
of environmental density (upper x axis). The horizontal dashed line 
indicates equal Mhi/M * > 1 and Mhi/M * < 1 galaxy frequency. 
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Fig. 23.— The frequency of HI gas-to-stellar mass ratios greater 
than one for ECO+A galaxies alone, plotted as a function of group 
halo mass for central (thick lines) and satellite (thin lines) galaxies 
separately. The horizontal dashed line indicates equal Mhi/M* > 
1 and Mhi/M* < 1 galaxy frequency. Confused sources in the 
ECO+A sample have been omitted. Since ECO contains rela¬ 
tively few satellite galaxies in M^aio < 10 11,5 M© environments, 
we refrain from plotting satellite frequencies in this regime. Since 
ECO contains relatively few centrals above Mh a i 0 > 10 13-5 Mq, 
we likewise refrain from plotting central frequencies in this regime. 
Gas-dominated galaxy fraction increases significantly for low group 
halo mass central galaxies, particularly at M^aio < 10 11,5 Mq. 

This offset may imply that for Mb ary < 10 10 Mq galaxies, 
disks are either destroyed less frequently or regenerated 
more frequently than for higher mass galaxies. 


5.1.2. Morphology-Environment Relations and Disk 
Regrowth 

We next consider the specific question of whether or 
not morphology-environment relations operate in a man¬ 


ner consistent with the presence of disk regrowth. If mor¬ 
phological transformation operates primarily towards the 
destruction of disks in certain regimes but towards both 
the destruction and regrowth of disks in others, we would 
expect the balance of galaxy morphological types to dif¬ 
fer in these regimes. In a scenario where large-scale gas 
accretion, whether arriving cold or hot, can fuel disk re¬ 
growth, the significance of such accretion is typically the¬ 


galaxy resides fe.g. 

Birnboim & Dekel|2003 

Keres et al. 

2005 Nelson et al. 

2013 

). Thus, in such a scenario, the 

balance ot galaxy morp 

nological types might naturally 


be expected to shift as a function of group halo mass. 
As previously mentioned, in the traditional morphology- 
environment relation, P(M|E), we observe a changing 
balance between early and late types as a function of 
group halo mass, in the sense that late types become 
more prevalent with decreasing halo mass. This trend 
is the sense in which the early/late type balance would 
be expected to vary if disk regrowth were to preferen¬ 
tially occur at low group halo mass, however it is also 
certainly the sense in which one would expect the rela¬ 
tion to vary if disks are typically destroyed/quenched at 
high but not low group halo mass. Thus, while consis¬ 
tent with a halo-mass dependent disk regrowth scenario, 
the observed traditional morphology-environment rela¬ 
tion does not clearly constrain its existence. We also 
note that the balance between early- and late-type fre¬ 
quencies we observe in the ECO sample varies relatively 
smoothly as a function of group halo mass, as has often 
been observed by other authors. This smooth variation 
could indicate a lack of sharp transitions in the onset 
of morphological transformation processes at particular 
mass scales, but alternatively it could imply that the tra¬ 
ditional relation, in lumping all early and all late types 
together, washes out possible sharper trends that may 
occur for subpopulations of galaxies. 

Considering this question from the perspective of an al¬ 
ternative formulation of the morphology-environment re¬ 
lation, P(E|M), another expectation of the disk regrowth 
model emerges. If disk regrowth were to proceed from 
blue early to blue late types in a particular regime, then 
the typical environments of blue early and blue late types 
in that regime should be similar as these galaxies would 
represent snapshots of pre- and post-transition states. 
Such behavior was hinted at in the observation of a pos¬ 
sible “inverse morphology-density relation” at low stel¬ 
lar masses by KGB. As illustrated in Figs. |20| and [2T| we 
find that blue early and blue late type galaxies in the low 
baryonic mass regime below ~ 10 1o Mq inhabit environ¬ 
ments with similar typical group halo masses at constant 
baryonic mass and with typical environmental densities 
of blue early types similar or lower than those of blue 
late types. The P(E|M) formulation of the morphology- 
environment relation then appears to be consistent with 
the scenario of disk regrowth in the low baryonic mass 
regime. If gas accretion adds to both galaxy baryonic 
and overall group halo mass during the regrowth pro¬ 
cess, individual galaxies could move along the Mhaio vs. 
Mbary relation, but relatively significant changes in typ¬ 
ical population properties would be necessary to make 
the blue early and blue late type populations distinct 
in this space given the scatter within each population. 
The trend towards lower number density environments 
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Fig. 24.— Characteristic distribu tion of HI gas-to-stellar mass ratio as a func tion of group halo mass for different galaxy types. HI gas 
masses are derived as described in §3.7| Symbols and lines are analogous to Fig. |20| 


for blue early types compared to blue late types possibly 
points to a post-merger status for blue early types. In the 
spirit of the alternative formulation of the morphology- 
environment relation, consideration of the typical envi¬ 
ronments of various subclasses of galaxies leads to fur¬ 
ther insights regarding environmental thresholds in disk 
growth as discussed in the next section. 

5.2. The Regime of Extreme Gas Richness and Recent- 
Disk Growth 

We have presented multiple results that add up to the 
impression that extremely gas rich galaxies and those po¬ 
tentially regrowing disks are preferentially found in the 
low group halo mass ( M ha i a < 10 115 M 0 ) and low galaxy 
mass (.Mbary or M* < 10 1o Mq, near the “gas-richness 
threshold” stellar mass at M* ~ 10 9 7 Mq) regime. We 
note that this regime is approximately defined, as in some 
cases we observe a continuum of galaxy properties in 
our sample, whereas in others, typically involving cen¬ 
tral galaxies, we see more abrupt transitions in properties 


between well-defined regimes. In Fig. 28 we summarize 
the galaxy mass and group halo mass distributions of 
the relevant galaxy populations. In ECO, blue-sequence 
early types and early types with Mm/M* > 1 occur 
most commonly below a group halo mass of ~ 10 n ' 5 Mg. 
This group halo mass regime is also where extreme gas- 
to-stellar mass ratios commonly emerge in our sample, 
with a relatively sharp uptick as a function of decreasing 
group halo mass as can be seen in the lower right panel 
of Fig. [28] Likewise, early-type UV-B galaxies are most 
common in the same low group halo mass and low galaxy 
mass regimes, albeit with a broader distribution extend¬ 
ing to somewhat higher masses as well. Early types with 
UV-B disks are more often gas rich than early types with¬ 
out UV-B disks, and most low-mass blue-sequence early 
types in the ECO+G subsample host UV-B disks. 

A possible explanation for this constellation of results 
is that this low galaxy mass and low group halo mass 
region of parameter space represents a preferred regime 
where gas is abundantly available to galaxies, fueling star 
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Fig. 25. — Group halo mass distribution for blue early type galax¬ 
ies in the ECO sample with different levels of HI gas. Correspond¬ 
ing greyscale dashed lines indicate the distribution for group cen¬ 
tral galaxies of each type. HI masses are derived as described in 
§3.7| Gas-dominated blue early types and blue early types with 
more moderate gas-to-stellar mass ratios primarily inhabit envi¬ 
ronments with group halo masses below ~ 1 0 12 Mq . with ~53% 
of gas-dominated blue early types in M^alo < 10 11 ' 5 Mg environ¬ 
ments. 


formation that allows many early types to live on the blue 
sequence, develop UV-B disks, and potentially even re¬ 
grow larger disk structures. The existence of such a gas 
and star formation rich regime could be a symptom of a 
large-scale cosmological accretion process that is partic¬ 
ularly efficient at supplying gas into galaxies at low mass 
scales. One such theorized process is “cold-mode” gas ac¬ 
cretion, thought to preferentially act at group halo mass 
scales below this ~ 10 n - 5 M© mass scale at z (e.g., 
the 10 11 ' 3 — lO 115 Af 0 scale of Keres et al. 12005 Keres 

iha 


et al.]2009). This model has recently been challenged by 


simulation results using the AREPO code, which imply 
that cold-mode accretion is not as significant for galaxies 
residing in low mass halos as previously thought, how¬ 
ever these results suggest that accretion of heated gas is 
a more significant contributor in this regime, causing the 
level of gas accretio n into low mass halos to rema in high 
([Nelson et al.|2013 ). As seen in the simulations of |Zehavi 


et al. 1 2(J12|) , from £ ~ 1 to the present the bulk of the 


stellar mass growth in low-mass halos (Mhaio 10 12 Mq 
is still due to star formation, while in high-mass halos, 
such growth is mainly due to mergers. 

Within the galaxy baryonic mass range we consider, we 
find that blue early-type and UV-B disk host galaxies are 
typically centrals in the Mhaio < 10 n - 5 M Q regime (see 
Figs. |Tf|) and 27). We also find that the strong uptick in 
gas-dominated galaxy fraction in this regime is primarily 
a central galaxy phenomenon, which may point towards 
accretion fue ling of central galaxies in these environments 
(see Fig. [23]). However, we also note that in the lowest 
group halo mass environments we probe there are rel¬ 
atively few satellite galaxies within our baryonic mass 
range, and therefore it is plausible that gas-dominated, 
disk-growing satellites with lower masses could be com¬ 
mon in such environments as well. From simulations, 
typical z ~0 gas accretion rates for satell ite galaxies m ay 
be lower than for central galaxies (e.g., Keres et al. 2009), 
but gas accretion may still pl ay an evolutionarily signif- 


icant role for satellites (e.g., Dekel & Birnboim 2006 


Simha et al.|2009). 


6. CONCLUSIONS 

In this work, we have considered two primary galaxy 
samples, the Environmental COntext (ECO) catalog, 
and the B-semester region of the REsolved Spectroscopy 
Of a Local VolumE (RESOLVE) survey. Both samples 
reach into the high-mass dwarf galaxy regime and span 
a variety of environments, with the larger ECO catalog 
sample including the greatest environmental diversity. 
Through comparison to the more complete RESOLVE-B 
catalog, we apply corrections for incompleteness effects in 
ECO, creating an approximately baryonic mass limited 
catalog down to 10 9 ’ 3 M©. In this analysis, we have em¬ 
ployed high-quality, custom-reprocessed optical, near-IR, 
and UV photometry along with morphological classifica¬ 
tions, atomic gas mass estimates, and multiple metrics 
of galaxy environment. 

Our key results are as follows. 

• We observe a traditional morphology-environment 
relation, P(M|E), similar to the expected form but 
with offset amplitudes between the low and high 
baryonic mass galaxy samples in the sense that late 
types are more common at low mass. 

• We find the form of the traditional morphology- 
environment relation to be consistent with the sce¬ 
nario that morphological transformation from early 
to late types (disk regrowth) could occur in a pre¬ 
ferred low group halo mass regime. However, this 
relation does not strongly constrain the existence of 
this scenario, as the relation can also be explained 
by the occurrence of disk destruction/quenching at 
high but not low group halo mass. 

• We consider an alternative form of the morphology- 
environment relation, P(E|M), which is instructive 
as a way of quantifying the typical environments of 
galaxies of various classes. This formulation leads 
to the observation that typical blue early-type and 
blue late-type galaxy group halo masses are simi¬ 
lar at constant baryonic mass, which is again con¬ 
sistent with expectations from the disk regrowth 
scenario. Likewise, the typical environmental den¬ 
sities of blue early types are similar or lower than 
those of blue late types at constant baryonic mass, 
potentially reflecting a post-merger state for blue 
early types. 

• The P(E|M) formulation of the morphology- 
environment relation also reveals that for central 
galaxies, there is no discernible relationship be¬ 
tween group halo mass and morphology at fixed 
galaxy mass: the typical halo masses for all early- 
type and all late-type centrals of a given mass 
are the same. However, satellite galaxies in dif¬ 
ferent color classes divide strongly in their typi¬ 
cal group halo masses at fixed galaxy mass, where 
red early and late types occupy higher group halo 
mass environments than both blue early and late 
types. These observations suggest the traditional 
morphology-environment relation is largely driven 
by a relationship between morphology and galaxy 
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Fig. 26.— Illustration of the typical galaxy and group halo mass distributions for ECO galaxies with HI gas-to-stellar mass ratios greater 
and less than one. HI gas masses are derived as described in §3.7| At fixed baryonic mass, high gas fraction satellites tend to occupy lower 
mass group halos than do low gas fraction satellites, implying group halo mass affects satellite gas content more strongly than galaxy mass. 
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Fig. 27.— Frequency and gas content of early-type UV-B disk 
hosts in the ECO+G sample. The magenta line indicates the fre¬ 
quency of UV-B disks in early types as a function of group halo 
mass, and the grey dashed line indicates UV-B disk frequency 
among early-type centrals alone. The inset shows the distribu¬ 
tion of HI gas-to-stellar mass ratios for UV-B and n on-U V-B early 
types. HI gas masses are derived as described in §3.7| A strong 
preference towards higher gas-to-stellar mass ratios is observed for 
early-type UV-B disk hosts. 

mass for centrals and by a relationship between 
color and environment for satellites. 

• We find that the low group halo mass regime be¬ 
low ~ 1 O 115 M 0 is associated with the emergence of 
blue-sequence early types, gas-dominated galaxies, 
and early-type UV-Bright disk hosts as common 
contributors to galaxy populations. These three 
sub-populations are closely linked in this regime, 
implying the low group halo mass regime is a pre¬ 
ferred regime for ongoing, significant disk growth. 

These results lend strong support to the idea that the¬ 
orized morphological transformation from early to late 
types can occur, particularly where galaxy and group 
halo masses are low. Thus, even if galaxy disks are de¬ 


stroyed through mergers and interactions in this regime, 
they may have significant opportunity to regrow, whereas 
galaxy disk regrowth does not appear likely in higher 
group halo mass environments. To investigate even more 
direct signatures of disk regrowth, we next turn to the 
examination of detailed early-type galaxy kinematics in 
the context of the RESOLVE survey, where the avail¬ 
ability of such kinematic information combined with the 
type of environmental information considered here cre¬ 
ates a unique opportunity for understanding the connec¬ 
tion between galaxy properties on small and large scales. 
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Fig. 28. — Illustration of the frequency (percentage) of various galaxy populations with the potential for disk growth as a function of 
galaxy and group halo mass. The top panels indicate the frequency of blue early- and late-type galaxies as a fraction of all galaxies with 
the same morphology and of all types. The lower panels show the frequency of gas-dominated and UV-B disk galaxies among all galaxy 
types and early types specifically. While each of these populations clearly displays a unique variation in frequency with both galaxy mass 
and environment, each of these “growing” populations occurs most commonly where galaxy masses and group halo masses are low. 
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